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Abstract. We review various theoretical approaches to the states of QCD matter 
out of quarks and gluons in extreme environments such as the high-temperature states 
at zero and finite baryon density and the dimensionally reduced state under an intense 
magnetic field. The topics at high temperature include the Polyakov loop and the 
't Hooft loop in the perturbative regime, the Polyakov loop behaviour and the phase 
transition in some of non-perturbative methods; the strong-coupling expansion, the 
large- iVc limit and the holographic QCD models. These analyses are extended to hot 
and dense matter with a finite baryon chemical potential. We point out that the 
difficulty in the finite-density problem has similarity to that under a strong magnetic 
field. We make a brief summary of results related to the topological contents probed by 
the magnetic field and the Chiral Magnetic Effect. We also address the close connection 
to the (1+1) dimensional system. 
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1. Introduction 

The inside of heavy nuclei is already a very interesting and peculiar environment. It is 
known that the density of nucleons, i.e. the baryon density, takes an almost constant 
value, po — 0.17nucleon/fm^, in the central region of nuclei independently of the atomic 
number A for large enough A. In terms of our daily units this normal nuclear density 
is as huge as ~ lO^^g/cm^. When two heavy nuclei (positively charged ions) collide 
at almost the speed of light, an enormous energy is crammed in a volume of size of 
heavy nuclei with the (transverse) radius ~ 1.2 A^^^ fm. In this way the relativistic 
heavy-ion collision experiment provides us with an ideal opportunity to examine the 
state of matter under extreme environments as have ever existed in the Universe. 

The collision energy is finally released into a form of the energy conveyed by 
produced particles. Measuring the momentum distribution of those particles, the initial 
energy density can be deduced, which is summed up as the Bjorken formula (though it 
looks slightly different from the original form [1]), 

^0 — — 2" ■ ~T~ ' y^) 
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where r = ^/t^ _ ^2 g^^^ V = \ ^^^[(^ + " ^)] are the proper time and the coordinate 
rapidity, respectively. Each produced particle has energy m_|_ = + rn^ with the 
transverse momentum p^. The number of particle is represented by A^. The energy 
density is then given by the energy of particle (mj,) diV divided by the initial volume 
TQ'Kr\dy, where tq denotes the initial time when particles are produced. 

Once the thermal equilibrium is reached and deconfined gluons and A^f-flavour 
quarks - the fundamental objects in Quantum Chromodynamics (QCD) - are assumed 
to be massless particles, the energy density e is translated into the temperature T 
through the Stefan-Boltzmann law; e = (vr^/SO) ri(T)T^ where n{T) = 16 (gluons) + 
lO.SA^f (quarks) is the effective number of physical degrees of freedom. The Monte- 
Carlo simulation of finite-T QCD discretized on the lattice (see references [21 [3] for 
recent advances) has identified the pseudo-critical temperature as Tc ~ (150-160) MeV. 
This value of Tc corresponds to the critical energy density, ec — (0.8-1.0) GeV/fm^, 
beyond which the state of matter should be composed of gluons and quarks, namely, 
a quark-gluon plasma (QGP) is realised. We note that the energy density at normal 
nuclear density po is ^nPo — 0.16 GeV/fm^ with the nucleon mass = 0.94 GeV 
used, that is about one fifth smaller than ec. 

In the facilities called the Alternating Gradient Synchrotron (AGS) at BNL and 
the Super Proton Synchrotron (SPS) at CERN the heavy-ion (Au-Au at AGS and Pb- 
Pb at SPS) experiments with a fixed target had been conducted since 1986. There, at 
the collision energy (per nucleon-nucleon) ^/s^ = 17.2 GeV at SPS, some theoretical 
studies |1] led to an estimate eo = (1.2-2.6) GeV/fm^ (at tq = 0.8 fm). This initial 
energy density exceeds the critical value, so that the state of matter created at SPS 
could be a QGP possibly. It has become evident in the Relativistic Heavy Ion Collider 
(RHIC) at BNL, which has been in operation since 2000, that the collision energy 
y/s^j^ = 200 GeV is high enough to form the QGP with the initial energy density 
Co = (12-20) GeV/fm^ (at tq = 0.6 fm). 

The experimental activities are still continued; there are two major directions as 
the future plan; one is the direction towards higher T at larger collision energies, while 
the other is the direction towards higher baryon density at smaller energies. The former 
is already ongoing since 2010 in the Large Hadron Collider (LHC) at CERN, where 
the experimental data from the Pb-Pb collision at ^/s^ = 2.76 TeV are significantly 
improving the quality of analysis. It is said that the QGP physics has been promoted 
from the "discovery stage" to the "precision science" studies. The latter direction, i.e. 
the extrapolation to higher baryon density regions and the experimental survey over the 
whole QCD phase diagram, is also underway. The beam-energy scan program at RHIC 
aims to give a detailed portrait of the phase structure of QCD matter, especially to 
locate a special point of the exact second-order phase transition called the QCD critical 
point [HI El |71 El [9] , which would serve as a landmark. The RHIC energy scan will be 
complemented by future experiments planned in the Facility for Antiproton and Ion 
Research (FAIR) at GSI, the Nuclotron-based Ion Collider Facility (NICA) at JINR 
and perhaps the Japan Proton Accelerator Research Complex (JPARC) at JAEA and 
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This review is a self-contained summary of some selected approaches in theory 
to the physical properties of QCD matter in extreme environments such as the high 
temperature in section |2] and the finite baryon density in section |3l It is, however, 
practically impossible to cover all the topics related to finite temperature/density QCD 
here. We shall specifically focus on the deconfinement physics and the dynamics of the 
order parameter called the Polyakov loop. As for the phase structure associated with 
chiral-symmetry breaking and restoration, the interested readers may consult my reviews 
[T0| [TT] and other reviews [121 [131 ttH [151 [IHl [13 [H] references therein. In particular 
my previous reviews [TOl [TT] include a pedagogical introduction to more generic physics 
of hot and dense QCD and subjects related to chiral-symmetry breaking, but not much 
about deconfinement physics. Hence, this present review is complementary to [lOl [TT] 
in respect to deconfinement physics and the Polyakov loop dynamics. 

The last half of this review is devoted to a new physics possibility in a strong 
magnetic field B produced by non-central collisions. The strength of this produced B 
surpasses the surface magnetic fields on the neutron star by orders of the magnitude. 
The presence of the magnetic background which is as strong as the QCD energy 
scale Aqcd ~ 200 MeV may well affect experimental observables considerably. A 
systematic procedure to resum relevant diagrams important at strong B has not been 
better established on the practical level than the finite temperature/density field theory. 
Therefore our discussions in section H] shall not be conclusive enough but they will aim to 
be comprehensive over various aspects of the phenomenon known as the Chiral Magnetic 
Effect. 

2. High-temperature State of QCD 

At sufficiently high temperature we can make use of the perturbative expansion in terms 
of the strong coupling constant, g, since the renormalisation of the UV divergences and 
the independence of the renormalisation point make g oi as = g^ I^t^ run as a function 
of the momentum scale [i\ at the one-loop order. 



set /io as the boson mass, Mz = 91.2 GeV and the world average at present is 
Q!s(M|) = 0.1184(7) [19]. It is natural from equation ([2|) to consider that the strong 
coupling constant gets smaller at higher temperature where the typical momentum scale 
among thermally interacting particles is characterised by the temperature T. One could 
therefore take n ccT but cannot precisely fix the proportionality coefficient because it is 
not clear which renormalisation condition is the most efficient to resum the higher-order 
diagrams. The conventional prescription is to take /i = 27rT and vary /i to check the 
stability of the physical results [201 (see also [2T] for a review). 



where /3o 
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Figure 1. (a) Graphical representation of the Polyakov loop in Euclidean space- 
time. The centre transformation multiplies on the Polyakov loop matrix, (b) Left: 
Polyakov loop correlator that is a counterpart of the closed Wilson loop. The correlator 
is centre invariant because zj. • = 1. Right: An opened- up figure whose temporal 
edges are contracted by the periodic boundary condition. 



While the perturbative QCD calculations are useful at high temperature (for the 
state-of-the-art calculations up to three- loop order, see |22l |23]), it is necessary to 
develop a non-perturbative method to go down towards T^. In this section we first 
look over the perturbative results, and then, we will proceed to several non-perturbative 
approaches. As we have stated, in this review, we mostly address the pure gluonic sector 
and the deconfinement order parameter. We postpone the discussions on dynamical 
quark effects to section [31 

The conventional choice of the order parameter for quark deconfinement at finite 
temperature is the (traced) Polyakov loop [211 125] . The Polyakov loop matrix and its 
traced quantity are denoted respectively as 



Pexp 



/3 



* = ^(trL), (3) 



in the imaginary-time formalism of the finite-temperature field theory. Figure [H (a) is 
a graphical representation on the manifold of §^ x M^. The Polyakov loop expectation 
value, $, can be interpreted as the partition function in the presence of a static-quark 
source. The logarithm of $ thus yields the single-quark free energy /q as in standard 
thermodynamics; /q = — Tln$. In the quark deconfined phase $ and /q take a finite 
value, whereas $ — )■ and /q — )■ cxd in the quark confined phase. 

This behaviour of the Polyakov loop is understood from Wilson's standard criterion 
of confinement ^26j. As seen in figure [H (b) the Wilson loop on x amounts to 
the correlator of the Polyakov loop L and the anti-Polyakov loop L^. In the confined 
phase the Wilson loop shows the area law (in the absence of dynamical quarks) and the 
deconfinement phase results in the perimeter law, which can be interpreted in terms of 
the Polyakov loop as 

Confined Phase 

{W{C)) ~ e""-^(^) ^ {iiL\r -> oo) trL(O)) -> =^ $ = (4) 

Deconfined Phase 
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{W{C)) ~ e-"wi(^) ^ (trL^(r ^ oo)trL(O)) ^ (const.) ^ $ ^ , (5) 

where S(C) = /3r is the area enclosed by C and is the perimeter (see figure [T](b)). 

This criterion is related to the global symmetry of QCD that makes the Polyakov 
loop expectation value vanishing. To manifest this symmetry let us turn to the lattice 
approximation to equation (j3]), i.e. 

Nr-l 

L = Y\ U4,{x,X4 = an^) . (6) 

714=0 

Here U4 = e'^^""^^ is the temporal link variable. (Note that in equation ([6]) the time 
product is from the left to the right with increasing time, while later time comes to 
the left in equation (jS]).) Then, under a general gauge transformation, the link variable 
changes as 

U^{x) V{x)U^{x)V^{x + afi) . (7) 

Therefore, the following transformation belongs to a subgroup of the gauge 
transformation, 

U^i^x, Xi = a{N^ - 1)) -> Zk U4,{x, X4 = a{Nr - 1)) (8) 

for all X with an A^c x Nc matrix, Zk = diag(e^'''''/^", e^'''''/^^ . . . , e^'"'''/^") where 
k = 0,1, . . . Nc — 1. (See figures [1] (a) and (b) for the graphical representation of 
the transformation.) Because this Z^^ group is a centre of the SU(iVc) gauge group, 
the global symmetry under the transformation (jH]) is called centre symmetry and the 
Polyakov loop changes accordingly as L ^ z^ L. Hence, the Polyakov loop expectation 
value is an order parameter for the spontaneous breaking of centre symmetry. 

In general the expectation value of operator can be determined so as to minimise 
the effective potential. Therefore, the effective action r[$] or the effective potential 
V[^] would suffice to give the information on whether the system is in the confined or 
deconfined phase. 

2.1. Perturhative approaches 

Because the perturbative calculations are valid at such high temperature that gluons 
and quarks interact weakly, the deconfined phase should be favoured in the perturbative 
regime and $ ~ 1 should be concluded. This anticipation was first confirmed by in 
[271 128] ill the SU(2) and SU(3) pure Yang-Mills theories. In the perturbative calculation 
it is more convenient to formulate the effective potential not in terms of $ directly but 
the phases of L instead. 

In this article we shall limit ourselves to the simple case of colour SU(2) only. The 
generalisation to colour SU(3) is straightforward. Then, with an appropriate choice of 
the basis in colour space with which is diagonal; 
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Figure 2. SU(2) Weiss potential (up to the normalisation) as a function of q (rescaled 
Ai, i.e. the phase of the Polyakov loop). The domain wall configuration interpolates 
between two minima. 



The Polyakov loop matrix and its trace can be expressed in terms of q as 

in the SU(2) case. Then, performing the one-loop integration on top of the 
background, one can find the following expression; 

viS[q] = 2VT J + ln[l - e-^^^''-^] + ln[l - e-^^-^^^"]} , (11) 

where the Polyakov loop enters as an imaginary colour chemical potential. After 
dropping the zero-point energy and carrying the momentum integration out, one can 
arrive finally at the perturbative effective potential (namely the Weiss potential), 

K?M = ^?Ldi(i-ta)^ (12) 

It is obvious from equation (ITT!) that v}q [q] is a periodic function of q with the period 
1, that is the reason why gmodi appears in equation (fT2!) . We remark that the Weiss 
potential is conveniently expressed by means of the Bernoulli polynomials. Figure [2] is 
a sketch of this periodic potential with the horizontal axis of q. In view of this periodic 
structure, one may well think of a tunnelling process that interpolates two minima. We 
will come to this point soon later. 

The two-loop calculations were first attempted in references [291 [30], which turned 
out incomplete because the one-loop correction to the Polyakov loop had been missing; 
in the computation of the Polyakov loop, the tree-level relation, (cos(7rg)) ~ cos(7r(g)), 
is insufficient for the two-loop potential. In this way, a simple prescription to treat 
the A4-background as if it were the Polyakov loop itself has a potential risk of pitfall 
in the non-perturbative regime. This correction was taken into the SU(2) calculation 
first [31J and the SU(A^c) calculation later [32]. Although the complete expression for 
the SU(A^c) case has a complicated combination of the Bernoulli polynomials, it is 
reduced to a simple form in the SU(2) case. The effective potential of the two- loop 
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order is composed from two pieces; [q\ from the direct contribution of the two-loop 

(1') 

integration and Vp [q] from the Polyakov-loop renormahsation at one-loop order; 

v!:^[q] = Vt\q]+Vi%]. (13) 
They are individually calculated and the SU(2) results read; 



2 2 2 

9'modl(l ~ Q'modl) — -Q'modl(l ~ Q'modl) 



?modl(l " Q'modl)^ — ^q'modl(l — fcodl) 



(14) 
(15) 



3/34 

Adding these two potentials up, we find that the latter terms inside of the parentheses 
cancel out and the final expression simplifies as 

Vff = ^^^?modl(l - ?modl) • (16) 

Therefore, interestingly enough, the correction of the two-loop order modifies only the 
overall coefficient of the Weiss potential and does not alter the functional form of the 
potential. We can confirm that the perturbative vacuum at g = is indeed a minimum 
of the potential and then $ = +1 is concluded from equation ( ITOl) . which is not ruined 
by the correction (fT6l) . 

From the Weiss potential, in the case when quarks are absent, the stable vacuum 
is degenerate at g = n with n being an integer, reflecting centre symmetry. These 
minima are, in fact, to be connected by the centre transformation. For general q = n 
the Polyakov loop takes a value of (—1)". This degeneracy between $ = +1 and $ = — 1 
would be broken by the presence of quarks. 

The quark contribution to the Polyakov loop potential is known up to the two-loop 
order [32]. We explain only the one-loop result here. The one-loop integration reads, 
apart from the zero-point energy, 

Crk[9] = -4iVf V^T I ^ {ln[l + e-^^^'-'^] + ln[l + e'^^-^^] } , (17) 
for massless A^f flavours, which eventually amounts to 



quark 3^4 1,2 2 / modi 



2 2 / modi 



It is worth mentioning that one can recover the above functional form immediately 
by replacing q — )■ g/2 + 1/2 in the Weiss potential f ll2l) . In this replacement an 
additional term 1/2 comes from the quantum statistics (boson or fermion) with which 
the exponential term changes the sign. Also, the argument is g/2 instead of q because 
quarks belong to the colour fundamental representation, while gluons are in the adjoint 
representation. We will take a closer look at quark effects when we discuss the model 
studies at finite density in section 13.31 

Since the period is doubled as compared to the pure gluonic case, the potential 
has no degeneracy between $ = ±1. Thus, $ = 1 (or g = 0) is more favoured than 
$ = — 1, which is a consequence of the explicit breaking of centre symmetry caused 
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Figure 3. Schematic figure of tlie 't Hooft loop. A singularity on the surface S twists 
the boundary condition for the Polyakov loop by a centre element. 



by quarks. In the presence of dynamical quarks, therefore, the state at $ = —1 is a 
metastable vacuum. It is shown [33], however, that the metastable state has physically 
unacceptable properties in thermodynamics. Different Polyakov-loop domains, whose 
interface is the Zjv^ domain wall [Ml |35], are meaningful only in Euclidean space-time 
(see [361 EZ] for detailed arguments against the physical interpretation of the Z domain 
wall). For the physical interpretation of the Polyakov loop in Minkowskian space-time 
the 't Hooft loop as we will discuss below is of special importance [38] . 

So far, we have considered only homogeneous configuration of the Polyakov loop 
background. One interesting application of the Weiss potential is the formation of the 
domain wall that is an inhomogeneous object in space. It is then convenient to introduce 
what is called the 't Hooft loop [39] (which measures the chromo-electric fiux) besides 
the Wilson loop (which measures the chromo-magnetic fiux). We denote the 't Hooft 
loop along the contour C as V{C) and then it should satisfy the operator relation; 

V\C)W{C')V{C) = e^i-Lkico/N.y^^^^;^ (^g^ 

with a centre-element coefiicient in the right-hand side and Gauss' link number Lk{C, C) 
of the two contours C and C . Mathematically Lk{C, C) can be expressed as [IQ] 



Lk{C, C') = ^f dx, I dy, e,,, , (20) 

J c J I *^ 2/ 1 

which shall be easily understood from Ampere's law with the magnetic field given by 

Biot-Savart's law in electromagnetism. From the fact that the 't Hooft loop counts the 

chromo-electric fiux, with a certain choice of the colour direction (which can be taken 

as the 3-rd direction in the SU(2) case without loss of generality), the explicit form of 

the 't Hooft loop could be given as [38] 

'27ri 



V{C) = exp 



d'S'Ef 



(21) 

9 J 

where S is the two-dimensional sheet enclosed by C, i.e. (9S = C. In the computation 
of the 't Hooft loop expectation value, the insertion of this operator to the functional 
integration induces a delta-function singularity (Dirac surface) on S, which makes a 
twist on the boundary condition for the Polyakov loop by a centre element (that is 
— 1 for the SU(2) group). We illustrate a schematic picture in figure [31 
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With the choice of C in the whole x-y plane, the surface E spans over the x-y plane 
at a certain (zq, ^o)- The expectation value of the 't Hooft loop is then expressed as 



{V{C)) = 1^ , (22) 

^pbc 

where Zpbc is the partition function with the periodic boundary condition in the z- 
direction, while Zthc has a twisted boundary condition for the Polyakov loop. 

Hence, to compute the expectation value of the 't Hooft loop, an effective action in 
terms of q is necessary including the derivative terms, that is given by 



(23) 



at one-loop order. This potential term is nothing but the Weiss potential and the 
derivative term is from the tree-level action. We note that the quantum corrections to the 
derivative term have been evaluated beyond the derivative expansion and a possibility 
of spatially inhomogeneous configuration of the Polyakov loop has been suggested |42] . 

In the presence of the 't Hooft loop, the twisted boundary condition is that g — )• 
at 2; — i- —00 and g — i- 1 at 2; — > +00. The 't Hooft loop is a creation operator of centre- 
domain interface, therefore. The classical solution associated with the above effective 
action is 

'^^^'^ = TT 1 ^ /^T^ ' ^^^^ 

1 + exp[-A/2/35(TzJ 

which satisfies the boundary condition; q{z — —00) = and q{z — > +00) = 1. Then, 
the effective action takes a finite value that is [34], 



47r2 

Tefffecl^:)] = ^-^=— L^Lj, = atL^Ly , (25) 

From this with equation (!22|) . the expectation value of the 't Hooft loop at the one- loop 
order is, 

{V{C)) = exp[-atS(C)] (26) 

with E(C) = LxLy, which shows the area law in the deconfinement phase. We see that 
the 't Hooft loop has behaviour opposite to the Wilson loop, and the 't Hooft loop plays 
the role of disorder parameter. Together with W{C) and V{C') the state of matter is 
characterised in more details as 

Confined Phase {W{C)) ~ e'^^^^^) , {V{C')) ~ e-<^^^'^ . (27) 

In the deconfinement phase or the Higgs phase, in contrast, the behaviour is opposite; 

Higgs Phase {W{C)) ~ e-<^'^^^ , {V{C')) ~ e'^^^^^'^ . (28) 

In the Higgs phase at high temperature, there is no massless particle. It should be 
noted that all gluons are massive due to the thermal screening mass. Here, we can also 
consider the third possibility; 

Partial Higgs Phase {W{C)) ~ e'^^^^^^ , {V{C)) ~ e'^'^^^^ , (29) 
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which represents the partial Higgs phase and confinement still remains. The last 
possibility that both the Wilson and the 't Hooft loops show the perimeter law is 
excluded from the operator algebra ( |T9i) [39]. 

Here we may understand that it is the spatial 't Hooft loop that exists in 
Minkowskian space-time and the real-time counterpart of the Zat^ interface is something 
that is created by the 't Hooft loop. The Polyakov loop in the real-time dynamics should 
be understood along this line |13] . 

We have, so far, discussed the behaviour of the Wilson loop or the Polyakov loop 
as the order parameter and the spatial 't Hooft loop as the disorder parameter. The 
spatial Wilson loop is also an interesting quantity. For completeness we shall give a brief 
description about the spatial Wilson loop. It always shows the area law regardless of 
the temperature. This can be understood in the 3-dimensional effective theory of QCD 
at high temperature [44j as a result of the dimensional reduction |45j . 

Integrating all "hard" modes out with non-zero Matsubara frequency at high 
temperature leaves a 3-dimensional effective theory of the "soft" length scales > {gT)~^. 
This effective theory is commonly referred to as Electrostatic QCD (EQCD), that is 
defined by the Lagrangian, 



^EQCD = j d3x|^trFi,F,,+tr(A^)' + m|trA2 + AE[tr(A2)2] + AEtrA 



.(30) 



The electrostatic field Aq{x) is static and adjoint scalar in colour space. The matching 
parameters in EQCD are calculated |46j and the three-loop level is still in progress 
[milH]. The leading-order results are; 

gl = g\T)T, ml = ^g\T)T\ Ae = , Ae = '\^J . (31) 

Because A^ix) is a heavily massive mode at high temperature, integrating out 
leads to an effective theory of the "ultrasoft" length scales > {g^T)~^, which is called 
Magnetostatic QCD (MQCD). This MQCD is defined by the Lagrangian, 



5mqcd = J d'x-ti F,jF,j . (32) 
This is a confining theory with the magnetic coupling constant, which is 

91a = 91 = g\T)T , (33) 

at the leading order. From the dimensional reason the string tension associated with 
this 3-dimensional effective theory (namely, the string tension measured by the spatial 
Wilson loop) is to be parametrised as 

(Ts = c^gl . (34) 

The determination of c requires full non-perturbative evaluation, and the Monte-Carlo 
simulation of the pure gluonic theory results in c ~ 0.553(1) [HI |50], which is also 
confirmed by later simulation [5T] . 

Although systematic resummation programs in EQCD and MQCD are on the track, 
the Polyakov loop effects and confinement physics are not incorporated in a satisfactory 
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Figure 4. Leading order contribution to the Polyakov loop effective action in the 
strong-coupling limit. 

manner (see reference [13] for some attempts and also [52] for another idea). It is quite 
difficult to investigate the nature of deconfinement phase transition in the perturbative 
approaches. 

2.2. Non-perturbative methods at work 

Theoretical researches on the confining properties near and below require non- 
perturbative extensions of the method. The lattice-QCD simulations are the most 
successful as long as the quark chemical potential is sufficiently smaller than the 
temperature. For recent developments in the lattice-QCD calculations there are a 
number of nice reviews (see reference [53] for example). In this review article we shall 
focus on some of analytical approaches. 

2.2.1. Strong- coupling expansion: The deconfinement phase transition can be 
formulated non-perturbatively in the limit of the strong coupling constant, — )■ 0, 
which was first elucidated in the Hamiltonian formalism in reference [24] . The same 
conclusion is readily obtained in the formalism of functional integration [Ml [55] . 

In the leading order of the plaquette expansion as sketched in figure HJ the effective 
action in terms of the Polyakov loop reads, 

= -e~'"'/^Y.^iL\xi)iiL{Xj) , (35) 

n.n 

which describes a hopping interaction between adjacent Polyakov loops. Here a is the 
lattice spacing. This action actually defines a spin-like theory of the Polyakov loop 
matrix; 

Z = j PLe-^p-'l^l . (36) 

Here T>L represents the functional integral with the group invariant (Haar) measure. 
The theoretical content of this matrix model itself is very intriguing [56]. In the same 
manner as the mean-field treatment (or the so-called molecular-field approximation) of 
spin systems, it is possible to formulate the spontaneous breaking of centre symmetry 
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and $ takes a finite value when the spin interaction becomes large at sufficiently high 
T [5711581 [591160]. 

Although it is much simpler than the molecular-field approximation, the tree- 
level potential is already useful to describe the deconfinement phase transition. In 
this prescription the traced Polyakov loop in the action (l35ll is simply replaced by the 
expectation value $, and an additional contribution comes from the Haar measure in 
the functional integration, i.e. 

Kfr[$] = -QVNl e-'^'^/^lxf - ln7WHaar[<f ] , (37) 
where the Haar measure for the SU(iVc) group is given by 

V ln[l - $$] , (for AT, = 2) 

V ln[l - 6$$ + 4(i>3 + $3) - 3($$)2] . (for A^^ = 3) 

It is important to note that these Haar measures favour the confining state at $ = 0. 
Moreover, the perturbative vacuum at $ = ±1 has an infinitely high barrier, which is 
cancelled by the longitudinal gluon loop in the perturbative calculation. Thus, the Haar 
measure could play an essential role in the realisation of confinement [611 ES] . Together 
with this Haar measure contribution and the spin interaction term, a phase transition 
takes place on the mean-field level and it is of second order for Nc = 2 and of first order 
for A'c = 3. Here, we distinguish the anti- Polyakov loop, $ = (ti L"^) /N^, from $; they 
are just identical at zero baryon density but a discrepancy between them arises from 
finite-density effects and has much to do with the sign problem. We shall return to this 
problem in the next section. 

The history of the investigations on chiral symmetry restoration in the strong- 
coupling expansion is as long as that of deconfinement physics, which is summarised in 
a review [63]. This section is devoted mainly to deconfinement physics and we will later 
look over the physics implications of chiral dynamics in section [31 

Inspired by the functional form from the strong-coupling analysis, one can adopt 
the following Ansatz to fit the pressure in the pure gluonic sector; 

V{<l>) = -^^<|$ + b{T) In 1 - 6$$ + A{¥ + $3) - 3($$)2 , (39) 

where a set of parameters, a(T) /T^ = 3.51 - 2A7t-^ + 15.2^-2 and b{T) /T^ = -1.75t-^ 
with t = T/Tc, can reproduce the lattice data well [HI]. In this parametrisation there 
are only three free variables because one of four is constrained by the Stefan-Boltzmann 
law. It is amazing that not only the Polyakov loop but also the pressure estimated from 

=0, P = V{^ = %), (40) 

simultaneously agree well with the lattice data [65] , as seen in the plots in figure [5l 

The agreement is impressive for only three fitting parameters, except for the 
Polyakov loop at very high temperature. The Polyakov loop from the lattice data 
exceeds unity there, which is caused by the renormalisation effect on the Polyakov loop 
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Figure 5. Left: Bands represent the lattice data for the Polyakov loop and the solid 
curve is the fitted result from equation (|40|) . Right: Pressure, the entropy density and 
the internal energy density from lattice data and the fitted result. The lattice data is 
taken from reference |65j . Similar figures are found in [66j . 

[671 168| [69| [65] . It is known that the U V divergence in the Polyakov loop is absorbed 
by the charge renormalisation as 

^'■^■^(T) = [Z(/)]^^$^^'^<=(/,A^,) . (41) 

The renormalisation constant is fixed at a reference temperature, that corresponds to 
the renormalisation condition. Once this is done in a certain scheme, a renormalisation 
constant at a different temperature is fixed by a different iV,-, which in turn leads to 
the renormalisation constant at a different coupling g (or lattice spacing). In this way 
the renormalised Polyakov loop is calculated for all temperatures through the iterative 
procedure. It is a non-trivial question how to incorporate the renormalisation effect in 
the Polyakov loop model |70] . 

2.2.2. Large-Nc QCD: The confinement- deconfinement transition is well-defined only 
in the pure Yang-Mills theories without quarks in the colour fundamental representation 
or in the limit of infinite quark mass. Otherwise, in the presence of dynamical quarks, 
centre symmetry is explicitly broken and the Polyakov loop always takes a finite value; 
$ 7^ 0. The Polyakov loop correlation function, in other words, does not decay 
exponentially at large distances due to pair creation of quark and anti-quark. 

Even in the presence of dynamical quarks, however, there is another limit in which 
the deconfinement transition is well-defined. That is, increasing the number of gluons 
instead of decreasing the number of quarks. In fact the quark contribution is more 
suppressed than gluons with large number of colours Nc ^ oo [711 [72], and eventually, 
in the limit of infinite colours a smooth crossover of deconfinement turns into a sharp 
phase transition. 

If we see the pressure of finite-T hadronic matter, on the one hand, it is of 0(1) 
in the counting. There are gluons of 0{N^), on the other hand, and the pressure of 
deconfined matter is of 0{N^). This is dominant over the quark contribution of 0{Nc). 
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Therefore, in the hmit of A^^c — !■ oo, the quark contribution becomes neghgible and the 
quenched approximation works correctly. The pressure then jumps from 0(1) to 0{N^) 
when the system goes through the phase transition from the hadronic to the deconfined 
phases. This means that the location of the phase transition has no ambiguity unlike 
the pseudo-critical temperature of crossover at A^c = 3. 

It is a subtle question what the order of the phase transition would be in the 
Nc oo limit. One might have thought that the phase transition should be of first order 
simply because of a big jump in the pressure. It is, however, possible to have a second (or 
higher) order phase transition with a pressure that is continuous but increasing rapidly. 
Thus we can naively think of two possibilities: (1) The phase transition (or crossover) 
is continuous for any A^^ and the ^ oo limit makes it of second order. (2) The phase 
transition at A'^c — !■ cxd is of first order and there is a critical number of colours at which 
the phase transition is of exact second order. 

The fact seems to be more complicated. In the large- A'c limit of the Polyakov loop 
matrix model another possibility has been suggested; the Gross- Witten point might be 
realised at A"c — )■ oo [73j. Then the effective potential is flat in the region < $ < 1/2 
and starts increasing for $ > 1/2. The Polyakov loop jumps from to 1/2 at the critical 
point. This is an unconventional point because no interface tension is needed for a jump 
between $ = and $ = 1/2, which would turn into a continuous transition immediately 
with an infinitesimal background. 

The large- A"c approach is useful as long as dropping dynamical quarks off is not 
critically harmful. It is an interesting theoretical challenge to apply the large-A'c 
argument for finite-density problems. 



2.2.3. Holographic model: The application of the AdS/CFT correspondence has 
become an important building block of hot and dense QCD physics. The idea is that the 
weak-coupling Type-IIB supergravity theory on AdSs x S5 is equivalent to the strong- 
coupling Af = 2 super Yang- Mills theory on the boundary of AdSs space. Especially the 
presence of the QGP is translated into a non-extremal black-hole solution (see references 
[711[75] for reviews). This technique is quite useful to examine non-perturbative aspects 
of strong-coupling gauge theories. The problem in the application to QCD physics 
is that QCD is neither conformal invariant nor supersymmetric. There are a number 
of theoretical attempts to design the black-hole solution so that it can mimic QCD 
thermodynamics [76l[771[78] (see also an extensive comparison with the lattice simulation 
[79] ) with a hope to estabhsh the AdS/QCD model. 

One of the simplest ways to introduce a mass scale is to use the following Ansatz 
for the five- dimensional background geometry (in Einstein frame); 

^-fiz)dt^ + dx^ + f-\z)dz^] (42) 



ds' = e^^ — 

with f{z) = 1 — (z/zii)'^, which describes a Schwarzschild-type black hole along the fifth 
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coordinate z. The corresponding Hawking temperature is 

T = — , (43) 

that is interpreted as the QGP temperature in the gauge theory side. The dilaton 
potential contains a dimensional parameter c that plays the role of QCD scale in this 
model. Such a model is usually referred to as the soft-wall model (usually defined 
in string frame) and is quite successful to give a semi- quantitative description of the 
Regge trajectory in the vector- meson channel [80] (see also reference [81] for a related 
approach). This type of approach is generally called the "bottom- up" model. 

Interestingly enough, the deconfinement phase transition is clearly identified in 
such a holographic setup. The QGP at high temperature is featured by the metric (H2|) . 
whereas confined matter at zero temperature is described by the metric with z\i oo 
(and thus T — )■ 0). Then, the five-dimensional actions associated with respective metrics 
determine which state is energetically favoured. In the soft-wall model the critical 
temperature Tc has been obtained as Tc ~ a/c/tt [H21 ES] where the Hawking-Page 
transition takes place. It is noteworthy that this phase transition of deconfinement is 
almost always of first order (see [77] for a possibility of exception). 

The Polyakov loop expectation value is calculated by the string world-sheet area 
that is minimised with the boundary along the Polyakov loop [HI]. It follows that 



$(T) = exp 



6|v^|Erfi(|)+l-e(^^/^)^ 



(44) 



where a is a normalisation constant and h = /2a' is a parameter in the Nambu-Goto 
string action. A choice of a = 0.10 and h = 0.72 fits the lattice data very well. 

As far as the Polyakov loop behaviour and bulk thermodynamic quantities are 
concerned, the results from the holographic approach are no better than the strong- 
coupling expansion as we have seen in equation fHU]) . It should be an advantage in the 
holographic model that some quantities that cannot be calculated on the lattice can be 
calculated easily such as the transport coefficients [851 ES] . Indeed, the computation of 
the shear viscosity in the strong-coupling expansion has not been successful particularly 
including the effect of the deconfinement transition 



3. Baryon-Rich State of QCD 

Historically speaking, the possibility of deconfined gluons and quarks was pointed out 
first for physics not at high temperature but at high density in the context of neutron 
star structure [88l[89]. It is, however, a non-trivial question whether the QCD running 
coupling constant really gets smaller at higher baryon density, as compared to the 
finite-T case in section [2l Because quark excitations are allowed only outside of the 
Fermi sphere if T is small, dynamical quarks must carry as large momentum as the 
quark chemical potential /iq. The important point is that the relevant scale in as{fi) in 
equation ([2]) is not the momentum of quarks but that of the exchanged gluon between 
quarks. It is still possible for fast-moving quarks to emit and absorb soft gluons, for 



QCD matter in extreme environments 



16 



200GeV 






lOGeV 


30GeV 










2Ge\\ 



200 400 600 800 1000 

Hb [McV] 



Quark-Gluon Plasma 

deconfmement 

X centre-sym 
light quarks 

O chiral-sym 




Color Superconductor 

deconfmement 
? centre-sym 
light CFL pions 
X chiral-sym 



T^fuclear Matter 



1^ 



B 



Figure 6. Left: Chemical freeze-out line. In the figure arbitrary five points are picked 
up to show the collision energy dependence of the chemical freeze-out point. Right: 
Minimal structure of the QCD phases in the ^b-T plane. The critical T is of order 
Aqcd and the critical /ie is near the nucleon mass which is of order TVcAqcd- 



which n could be small and as(/^) could be substantially large. Such soft processes are, 
however, screened immediately due to quark polarisation effects that induce a screening 
mass ~ (7/iq on gluons. In this way, in effect, one can regard sufficiently high-density 
matter of QCD as a weak-coupling system. This way of understanding is challenged 
recently by the large-A'^c approach to the QCD phase diagram as we will see later in 
section 13.2.21 

It is an ongoing experimental project to explore the state of QCD matter in a wide 
range of temperature and baryon density by varying the collision energy y/s^. From 
the phenomenological analysis using the thermal Statistical Model |90l |911 192] a set of 
the temperature T and the baryon chemical potential /ib (= 3/iq) at which the chemical 
composition of particle species is frozen has been extracted. These "chemical freeze-out 
points" are very well parametrised by [90] 

T{{fiB) = a - bfi'^ - cfi^ (45) 

with a = 0.166 ± 0.002 GeV, b = 0.139 ± 0.016 GeV^^ and c = 0.053 ± 0.021 GeV"l 
Also the baryon chemical potential at chemical freeze-out is parametrised as a function 
of the collision energy as 

'"<^' - TT7^ («' 

with d = 1.308 ± 0.028 GeV and e = 0.273 ± 0.008 GeV"\ from which one can convert 
equation f H5]) into Tf(ys^) easily. 

These parametrisations are useful to relate theoretical predictions (for example, 
fluctuations of the charged particles, the baryon number etc; see reference [93]) to 
experimental observable as a function of not only /ib but also y/s^ that is under 
experimental control. 

It is a natural anticipation that this chemical freeze-out line, which is plotted in 
the left of figure |6l is related to the colour deconfinement phenomenon that causes 
a rapid increase in the particle number density, and thus the multiparticle scattering 



QCD matter in extreme environments 



17 



rate [94] . We note that chiral symmetry restoration is not taken into consideration 
at all in the assumption of the thermal Statistical Model. In this way a part of the 
theoretically conjectured QCD phase diagram (the right of figure [6]) has been supported 
by experimental data, which should be complemented further by a line associated with 
colour superconductivity (see reference [95] for a modern review and references therein). 
We should emphasise that figure [6] shows only the minimal structure of possible QCD 
phases. There are many other possibilities on top of this minimal topology, some of 
which are summarised in my previous review [TT] . 

Progresses in the experimental activities towards QCD matter at higher baryon 
density are strongly needed by theorists. Unlike the case at high T and small /xb, 
theoretical works have not been successful in making any robust prediction on the 
baryon-rich state of QCD matter. It is absolutely necessary to constrain proposed 
theoretical possibilities from the experimental point of view. 



3.1. Perturhative approaches and problems 

Let us first consider an extension of the Weiss potential ( IT2l) to the finite-density case. 
We can perform the perturbative integration as in the standard procedure with a quark 
chemical potential /Xq introduced [96l |971 [98]. It seems at a first glance that such 
extension is straightforward, but it is not so simple once the Polyakov loop background 
is involved. 

The loop corrections to the Weiss potential have been evaluated with dynamical 
(massless) quarks. The one-loop contribution from massless quarks is modified by /iq 
from equation (ITTIl into 

KJuLM = -2iVfVT j 5^{ln[l + e-PP+P^^^+^^-<i^] + ln[l + g-Z^P-Z^Mq-i-^,] | ^ (47) 

which is generalised from SU(2) to SU(A^c) with j running from 1 to (and qj = 0). 
The first and second logarithms represent the particle and the anti-particle excitations, 
respectively. Then, the result after the momentum integration is obtained immediately 
[99] by the replacement of 

9. ^ - (48) 
in equation (ITS!) . 

This final result is simple but astonishing. Unless the gauge group is SU(2), 
the effective potential generally takes a complex value. How can one determine the 
energetically favourite value of qj from such a complex potential? One might have 
thought that qj should minimise the real part of the potential. Although such a 
working hypothesis may give a practical prescription, this cannot be justified from the 
first-principle approach. This complex potential for the Polyakov loop is one clear 
manifestation of the notorious sign problem (see reference [100] for an introductory 
review) . 
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The sign problem hinders the lattice-QCD simulation at finite density. Contrary to 
what is believed, the sign problem is actually a quite generic problem of the importance 
sampling not only in the lattice-QCD simulation but in the mean-field approximation 
also |10H 1102] . One should notice that the mean-field approximation is based on the 
importance sampling; the mean-field variables are chosen to be a "configuration" that 
maximises the weight ~ Q-'^c,ft[q]_ When the potential is complex, therefore, the mean- 
field approximation breaks down. 

Occasionally, in some analytical studies, it is overemphasised that the method be 
sign-problem free. Such a statement must be misleading as long as the method relies 
on the mean-field approximation for the treatment of the gauge-field part such as the 
Polyakov loop dynamics. 

3.2. Non-perturbative methods in progress 

One of the most urgent challenges in theory is to outline the global structure of the QCD 
phase diagram and fill in the quantitative details on figure El Non-perturbative methods 
are indispensable to access the information in the vicinity of phase transition regions. 
There are significant progresses recently in the strong-coupling expansion, the large- 
Nc QCD and the effective models to shed light on the phase diagram. Unfortunately, 
not much about the phase diagram can be said from the holographic QCD models, 
though there are many interesting attempts on each state of QCD matter, particularly 
by means of the AdS-Reissner- Nordstrom black hole |103l 11041 11051 1106] . where the 
charged black hole in AdS space is identified as the finite-T and finite-/iq plasma. There 
are some more investigations on quark-matter and nuclear-matter properties using the 
holographic approach [107[ 1108] and also on the phase transition in the Sakai-Sugimoto 
model [1091 IllOl llllj . In this article, not only the holographic approaches but also the 
finite-density studies by means of the orbifold equivalence |112[ 1113] and the orientifold 
equivalence |114l 1115] are beyond our current scope. 

3.2.1. Strong- coupling expansion and the matrix model: The Polyakov loop matrix 
model emerges as a result of the strong coupling expansion and this model provides us 
with an ideal setup to think of the sign problem. In the leading order of the hopping 
parameter expansion [116] , in the presence of heavy quarks, the quark- loop contribution 
or the Dirac determinant amounts to 



with a coefficient h which is small for large quark mass and eventually A^quark(/^q) 
1 as rriq — >■ oo (quenched limit). Because the Polyakov loop changes non-trivially under a 
centre transformation, the above quark contribution breaks centre symmetry explicitly. 
Besides, this action becomes complex when a finite /iq is turned on since trL takes a 
complex value in general. 




(49) 



X 
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The expression (H^ is very useful to grasp the nature of the sign problem. One can 
also use equation fH9l) the other way around to deduce special situations where the sign 
problem weakens. We shall enumerate some of widely acknowledged examples here: 

• Two-colour QCD [HIl [IHl [119] — If the colour gauge group is SU(2), the 
Polyakov loop is always real; trL = trL^ = 2cos(7rg). Equation ( l49l) is thus 
real. Taking a real value is not sufficient for feasibility of the importance sampling 
because the Dirac determinant could be real but negative. There must be an even 
number of degenerate quark species in order to guarantee the semi-positivity of the 
Dirac determinant. 

• Isospin Chemical Potential [120(1121] — We see that A^quark(— /Uq) = -^quark(Aiq)* 
from equation (H9l) . meaning that the whole Dirac determinant is positive semi- 
definite if there are two degenerate quarks that have a chemical potential opposite 
to each other. For example, = /ii for w-quarks and Hd = —fJ'i for (i-quarks with 
rriu = m^. It is then easy to confirm that the Dirac determinant satisfies, 

quark (/i«)A^ quark (/Ud) = | quark (/"l) P > . (50) 

The isospin chemical potential causes no sign problem, therefore, and the Monte- 
Carlo simulations are feasible. The absence of the sign problem for the chiral 
chemical potential /is is also a variant of this category |122] . for which some lattice 
simulations are successful |123[ 1124] . 

• Imaginary Chemical Potential |125[ I126[ 1127] — The sign problem originates from 
the imbalance between the quark and anti-quark propagation in equation f H9|) . 
which can be made balanced by replacing /iq by a pure-imaginary quantity i/iq. 
Then A^quark(Aiq) IS obviously real. In this case, unlike two-colour QCD, there 
needs not be an even number of degenerate quark species because e^'^'^'^ is bounded. 
(The situation about the positivity of the Dirac determinant is rather similar to 
the zero-density case.) One may have thought that the Dirac determinant is then 
a periodic function of /tq with a period 27rT. This is not correct. In fact the phase 
factor e'^^"! can be partially cancelled by the centre transformation and the rest 
takes a value from 1 to e^™/^'^. This means that the genuine period is 2-kT/Nc 
instead of 27r (i.e. Roberge- Weiss periodicity [128]). 

Even though the sign problem is not washed away, the mean-field approximation 
works anyway (see section l3T3] for details). In the strong-coupling expansion the gauge 
action is dropped and the quark sector is dominant. Then, it is not the deconfinement 
transition but the chiral phase transition that defines the phase diagram. The phase 
boundary of chiral restoration obtained by means of the staggered formalism of chiral 
fermions [631 11291 1130] is in qualitative agreement with figure [6l which has been also 
confirmed by the numerical simulation [131] . 

The effect of the Polyakov loop dynamics as formulated by the action fl35|) has 
been taken into account too [STJ [5H1 EHl [60]. Such a treatment on the lattice can be 
easily translated into the continuum language, which has led to the so-called PNJL-type 
models as we discuss soon below. 
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Figure 7. The phase diagram of large- iVc QCD ,132; . The hadronic phase confines 
coloured excitations and only glueballs and mesons exist. The corresponding system 
pressure is of 0{N^). Above the deconfinement temperature of order of Aqcd, thermal 
excitations are dominated by gluons and the associated pressure is of 0{N^). The 
quarkyonic phase is characterised by the pressure of 0{Nc) at greater than the 
lightest baryon mass so that the baryon number density is non- vanishing. 

3.2.2. Large-Nc QCD: Recently an interesting possibility about a new structure on 
the QCD phase diagram has been suggested from analytic deliberations on the large- A^c 
limit of QCD at finite T and /xb [132] . The phase diagram takes a simple structure as 
sketched in figure[7]with three regions separated by straight first-order phase boundaries. 

In the large- Ac limit quark loops are suppressed by 1/Ac as compared to gluon 
loops, so /iB does not affect the deconfinement transition line which is predominantly 
determined by A^ — 1 gluons. Hence the deconfinement transition makes a straight line 
parallel to the /^b axis. A free quark gas would give a pressure of ~ A^/iq ~ N~^fi^, and 
thus the deconfinement transition line would be deformed for /iB as large as 0{Nc^^) that 
surpasses the gluon degrees of freedom. Before this is reached, there appears another 
type of transition at /iB — where M^r is the lightest baryon excitation energy (baryon 
mass minus binding energy). The baryon number density becomes non- vanishing then 
and the system pressure jumps from O(A^) in the hadronic phase to 0(Ac) in the new 
phase which is called Quarkyonic Matter. Before the deconfinement phase transition 
takes place with T ~ Aqcd, glueballs and mesons cannot affect this threshold for the 
baryon number and so the threshold located at /iB — M^y makes a straight line parallel 
to the T axis, which results in the phase structure presented in figure [71 

The reason why the right-bottom region of figure [7] is identified as the Quarkyonic 
Matter is the following. As explained before, gluons and quarks are all confined below 
the deconfinement transition line, and so the Quarkyonic Matter region resides in the 
confined regime. Therefore the physical degrees of freedom there should be baryons 
rather than quarks, and one can prove that the dense baryonic system at large Ac 
indeed gives the pressure of ~ 0(Ac) whose major contribution comes from not the 
Fermi energy but the baryon-baryon interaction energy |132j . On the other hand, a gas 
of quarks naturally yields a pressure of 0(Ac) because of the presence of Ac quarks. Such 
coincidence in the Ac counting implies that this bottom-right state would confine quarks 
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and nevertheless feel quarks somehow. In other words a dual interpretation is possible; 
matter of strongly-interacting baryons and simultaneously that of weakly-interacting 
quarks, which motivated the name, Quarkyonic Matter = (Quark + Baryonic) Matter. 

Such an interpretation may sound peculiar but there is a reasonable way to reconcile 
two interpretations in terms of baryons and quarks. That is, particles sitting deeply 
inside of the Fermi sea cannot be a part of excitation spectra, and so they could be 
quarks even though the system is in the confined phase. Therefore, the Fermi sphere 
consists of both baryons and quarks; baryons in the momentum layer ~ Aqcd near 
the Fermi surface (but, because of the confining interaction, there is no sharp Fermi 
surface in reality) and quarks inside of the Fermi sphere which does not take part in the 
excitation but gives the pressure of ~ 0{N^. 

This is actually what should be expected in the \sige-Nc limit. As we have discussed 
in the beginning of this section, a large /iq does not guarantee the smallness of as(/Xq), 
but the screening effects due to quark polarisation make a^^ii) effectively small enough to 
realise the perturbative regime at high density. Because the quark loops are suppressed 
in the counting, the screening effects would diminish and soft-gluon exchange would 
become important for large A^^c. 

The soft-gluon exchange would lead to confinement for excitations on top of the 
Fermi surface, and furthermore, to an interesting consequence for the chiral-symmetry 
breaking mechanism |133l 1134] . Usually the homogeneous chiral condensate (gg) or the 
condensation of the sigma meson at rest is attributed to the spontaneous breaking of 
chiral symmetry. Because quarks and quark-holes must sit near the Fermi surface, if 
the net momentum of gg-system is zero on the one hand, the gluon exchanged between 
a quark and a quark- hole on the Fermi surface must carry as large momentum as 2yUq. 
On the other hand, this momentum can be absorbed in the net momentum p ~ 2/iq of 
gg-system and then the possible soft-gluon exchange favours an inhomogeneous chiral 
condensate, namely, the chiral spiral state. 

It is a subtle question whether the chiral spiral structure emerges and Quarkyonic 
Matter exists in the real world at A'c = 3. This depends on the competition of the 
strength of confining force and the screening due to quark polarisation. For the sake 
of such quantitative clarification the effective model study was expected to hint the 
relevance or irrelevance of the large- A'c limit to the real world. It has turned out, 
however, that one should inevitably go beyond the mean-field approximation and it is 
still a very difficult question how to resum the Polyakov loop fluctuations. 

3. 3. Extrapolation from effective models 

The Polyakov loop behaviour in the pure gluonic sector is described nicely by the 
simple parametrisation ( l39i) . The quark loop on top of the Polyakov loop background 
in equation (H7I) is expressed in a gauge invariant way, i.e. 

^quL = -mVT j ^|trln[l + Le'^(^-^-)] + trln[l + L'^ e-^^'+^-^]^ , (51) 
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where the quark mass is included in the energy dispersion relation e = _|_ ]\^2 -^j^Ji 
the dynamically generated mass Mq which could be considerably larger than /iq. This 
form of the coupling is simple but has rich contents. If chiral symmetry is badly broken 
by large Mq, the exponential terms are small and thus the Polyakov loop coupling 
to the chiral sector diminishes, which is in favour of confinement. If the Polyakov- loop 
expectation value is small, on the other hand, the thermal excitation of quarks is severely 
screened and thus chiral restoration is delayed. 

Interestingly enough, the colour trace in equation ( l5Ti) is explicitly taken to be a 
form of 

y^l^. = -mVT j ^|ln[l + trLe-^(^-^'')+trLte-2''(^-'^^) + e-3^(^-'^'') 

+ In 1 + tr e-^('+''^) + tr L e-2^(^+^<') + | . (52) 

The further missing piece in the dynamics is the spontaneous breaking of chiral 
symmetry now that the deconfinement and coupling parts are formulated as explained 
above. The Polyakov-loop coupled Nambu-Jona-Lasinio model (PNJL model) |135tll36] 
utilises the NJL model as a dynamical theory to describe the chiral condensate {qq). 
Replacing the NJL model by another successful chiral model, the Quark-Meson (QM) 
model (that is a variant of the linear sigma model), one can define the PQM model 
|137] . In the simple mean- field approximation the PNJL model is more convenient than 
the PQM model because the linear sigma model suffers from artificial first-order phase 
transition |138] . This model artifact is cured by pion loop effects implemented by the 
renormalisation group (RG) improvement. For this purpose of RG studies the PQM 
setup has an advantage over the PNJL model |139[ 1140] . 

So far the RG improvement on the meson fluctuations has been well investigated, 
while it is not known how to include the Polyakov loop fluctuations systematically. 
Once the Polyakov loop potential is given, there is no way to improve it, and a more 
fundamental starting point is necessary. One possibility is to take the Polyakov loop 
matrix model that can in principle encompass soft fluctuations of the Polyakov loop. 
Since this matrix model is defined on the lattice, however, it is technically difficult to 
accomplish the RG analysis. Another possibility is to postulate a derivative term in the 
Polyakov loop action, i.e. setup of the so-called Polyakov loop model |141] and to put 
it in the RG equation |142j . 

Here, let us make a brief remark on the sign problem in the PNJL-type models. 
It may seem to be sign-problem free in equation ( 152|) but it is not so. In the simple 
mean-field analysis equation (152|) leads to 

^quLf"^' <!] = -2NiVT j ^^|ln[l + 3$e-^(^^^'') + 3$ e-'^^^"'^^) + e'^^^^"'^'') 

+ In 1 + 3$ e-'^(^+'^<') + 3$ + Q-^Pi'+f'^) | . (53) 

Here, at finite /iq, it is necessary to treat $ and $ independently. Then, one can plot 
Kiuark[^' $] as a fuuctiou of $ and $ to recognise a funny shape. In particular one may 
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Figure 8. Entropy density and the baryon number density normalised by the 
Stefan-Boltzmann value in the mean-field approximation of the PNJL model. The 
entropy density, which is a T-derivative of the thermodynamic potential, should give 
information on gluon deconfinement, whereas the baryon number density, which is a 
/iB-derivative, should be sensitive to quark deconfinement. 

find that the solution of the gap equation is not stable in the direction of $ — $. This 
is how the sign problem remains unsolved in the mean-field model study and in fact 
identifying the saddle-point as the ground state leads to an approximation similar to 
the reweighting method to evade the sign problem |102j . 

Let us make a brief comment on the approximation: trL — )■ 3$ in the logarithm 
in equation f l53|) . Such a replacement is acceptable only when the Polyakov loop 
fluctuations are negligible as compared to its mean value $. To go beyond the mean- 
field level the group integration with respect to the Polyakov loop matrix L should be 
incorporated as done in [581 EHl EOl 1143] . Such a treatment of the group integration 
with respect to L is important to maintain gauge invariance; in the literature gauge- 
variant quantities such as the phases of the Polyakov loop are chosen as the mean-field 
variables to simplify calculations with colour-superconducting gaps, but they would lead 
to unphysical colour density. This problem is often neglected but must be resolved by 
the group integration as elucidated in |144j . 

Nevertheless, once an approximation is made with a prescription to neglect the 
sign problem, the model results are useful to figure the thermodynamic quantities out 
at finite T and /ie- Figure El shows some examples from the PNJL model in the simple 
mean-field approximation |145] . The left is the entropy density divided by the Stefan- 
Boltzmann value and the right is the baryon number density divided by the Stefan- 
Boltzmann value. Naturally the increase in the entropy density is to be interpreted as 
deconfinement. One may well conclude that the model results could have implied the 
realisation of Quarkyonic Matter at low T and high /xb where the entropy density stays 
small and the baryon density gets large which is characteristic to Quarkyonic Matter. 

However, the model study has missing diagrams, as depicted in figure [9l which is not 
included even in the RG improvement and is related to the Polyakov loop fluctuations. 
This missing contribution is critically important to clarify the QCD phase structure 
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Figure 9. Ring diagram of gluons with the quark pair creation/annihilation. The 
Polyakov loop potential must have coupling to the quark chemical potential /iq through 
this quark polarisation. If an external B is applied, gluons can feel B through this 
diagram too. 

including the possibility of Quarkyonic Matter. Without this diagram the Polyakov loop 
potential fl39|l has no explicit dependence on /iq, so that the deconfinement transition 
is almost insensitive to /iq even with the coupling effects through equation (|5T1) . This 
approximate /iq-independence is observed in the entropy behaviour too in the left of 
figure [HI Then, the phase diagram from the PNJL or PQM models turns out similar to 
the large- Ai'c conjecture in figured In fact, it is the polarisation effect in figure [9] that 
would make a difference from the large- A^^c limit. 

It is not straightforward to implement the polarisation effects properly in the 
effective model without losing simplicity. The model is to be appreciated as long as 
it is simple enough to deepen the intuition of physics understanding. There are only 
a few attempts to incorporate effects originating from the diagram in figure [9] into 
the PNJL-type models. By the hypothetical matching condition for the deconfinement 
and chiral transitions the /iq-dependence was introduced in reference [137] . Then, the 
deconfinement line comes along with the chiral restoration line on the phase diagram, 
but this statement is a consequence by construction of the model. In reference [146] 
the /iq-dependence in the Polyakov loop potential was determined by the matching 
condition to thermodynamics from the thermal Statistical Model, which has confirmed 
quantitative agreement with the prescription in reference |137j . 

This kind of polarisation effect is of increasing importance in the researches on finite- 
density QCD matter. Besides, as we see in the next section, QCD in strong magnetic 
fields is currently a hot subject and it requires detailed information on the polarisation 
effect. The reason for this is exactly the same as the finite-density case. Gluons do not 
feel the magnetic fields directly, but do see them through the quark polarisation as in 
figure [H 
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Figure 10. Left: Collision geometry seen from the above. Right: Profile of the 
produced magnetic field as a function of time. 



4. Strong Magnetic Field and Dimensional Reduction 

In the heavy-ion collision with a finite impact parameter (i.e. peripheral collision) a 
magnetic field is created by the positively charged ions moving at almost the speed of 
light. Let us evaluate how large magnetic field is expected in the collision at the RHIC 
energy in a classical manner. For simplicity we assume that the (positively charged) 
heavy ions are point charges |147j . The collision geometry is schematically modelled as 
in the left of figure [10 Then, from the Lienard-Wiechert potential, the magnetic fields 
at the origin reads 

. _ 6/3(1 -/32)ey _ Sy 

eB{x,t) - ■ ^^^^^2 ^ _ /32)(6/2)2]3/2 " [i + (t/t^ff/^ ' ^^^^ 

e5o = ^^sinh(F) = (47.6 MeVf Z sin\i{Y) , 

6 

° ~ 2sinh(y) ■ 

In the definition of Bq and to we use the beam rapidity Y instead of the velocity /3, 
which is related by /3 = tanh(y). Here, Bq is the maximum strength of the magnetic 
field and to gives a typical time scale of decaying field. In the case of Au-Au collision 
at the RHIC energy, these parameters are 



fs — 

Z = 79 , sinh(r) ~ cosh(r) = ~ 106.6 . (55) 

The point-charge approximation is valid when the collision is far peripheral. So, let us 
take 6 = 10 fm |147j . Then, this simple estimate leads to 

eSo ~ 1.9 X 10^ MeV^ = 3.2 X lO^Vuss , to ^ 0.05 fm/c . (56) 

This magnetic field strength is 10^ times larger than the surface magnetic field of the 
magnetar, and 10'' times larger than that of the ordinary neutron star. Although such a 
strong field is transient and decays with the time scale to? we note that the decay is not 
as steep as exponential damp but power-law suppression. At t/to ~ 2 ~ 0.1 fm/c, for 
example, the magnetic field diminishes to a tenth of Bq. We note that this time scale is 
of order of where Qs is the saturation scale at RHIC |1481I149[I150] . Although there 
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are a number of theoretical calculations of equilibrated QCD matter under strong B 
fields, any serious simulation of the Glasma |151j and the particle production in strong 
B have not been fully analysed. The real-time dynamics of the strong B effects needs 
more investigations. 

4-1. Topological properties probed by the magnetic field 

It is well-known that special gauge configurations with non-zero winding number play an 
important role in understanding of the vacuum structure in the strong interactions. The 
spontaneous breaking of chiral symmetry is attributed to the QCD instanton which is 
the origin of dynamical mass generation |152] . The confinement nature is also explained 
in terms of magnetic monopole condensation in a special class of the gauge choice (see 
[153] for a lecture note). 

There is no doubt about the existence of topological configurations in QCD physics, 
but it is quite challenging how to "see" such topological contents in real experiments. 
The Chiral Magnetic Effect (CME) is one of the promising candidates |122[ I154[ 1155] . 
Let us imagine the following situation specifically; the QCD vacuum accommodates one 
instanton that has a topological charge Qw and a magnetic field B that is as strong as 
the QCD energy scale Aqcd is applied on the instanton. 

Then, the axial anomaly relation (for the single- flavour case), 

d,j' = -£^ I d'xtTF,J^\ (57) 

implies that 

AN5 = N^it = oo) - Nr^it = -oo) = -2Qw . (58) 
The topological charge Qw is given by the temporal integration of the right-hand side of 
equation (1571) . This means that, if the system starts with the chirally neutral situation 
{N^lt = —00) = 0), a finite amount of chirality in the final state is generated by the 
topological charge during the time evolution. In the chiral limit the momentum and 
the spin are parallel to each other if the chirality is right-handed, while they are anti- 
parallel if the chirality is left-handed. The spin is aligned by strong B, which makes 
the momentum also aligned along the B direction, leading to a non-vanishing value of 
the total momentum if AA^5 7^ 0. In other words, since Dirac fermions are charged, an 
electric or baryonic current is produced for B ^ and AN^ ^ 0. 
Such an effect can be expressed simply as |154j 

Jv = -2QTy^ , (59) 

for large enough B, where Jy represents the vector current [ip'^^il)) integrated over the 
spatial volume. 

For arbitrary strength of B it is more appropriate to work in the grand canonical 
ensemble using the chiral chemical potential /^s instead of N^. One can then prove 
non-perturbatively that |122[ I156[ 1157] 

i. = ^B (60) 
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Figure 11. Chiral Magnetic Effect: In strong B field right-handed particles move in 
parallel with B and left-handed particles in anti-parallel to B. A non-zero net flow 
results from N^^O. 



holds for any B and /is. It is easy to confirm that equation ( l60i) is reduced to 
equation ( 159|1 in the strong B hmit using the anomaly relation (158|) . 
In fact this is one example of more generic relation, 

f = Ce^''"'''d,cpFp,. (61) 

Here (p is some field, and in the context of the CME, the strong 6'-angle is identified 
as <p in the above. Then, z/ = 0, p = x and a = y uniquely fix u = z from ei^'^c" _ By 
regarding OqO as /is |122j . one can readily retrieve equation ( 160|) apart from the overall 
constant C. One can also understand related effects from equation ( 16T]) . If 6 (or (p 
in the above) is spatially inhomogeneous, v = z (instead of z/ = 0) leads to /i = 
(instead of /i = 3), meaning that an electric charge or an electric dipole moment is 
induced |158] . In addition, one can even think of not the 6'-angle but the pion field tt 
as in equation (16T|) . and then the presence of time-dependent background, i.e. (9o7r°, 
results in an effect similar to the CME- induced current. It is actually argued in the 
Skyrmion model that an additional electric charge is induced in the baryon content 
(that is constructed as a profile of 7r"(a;)) under a strong B field |159j . 

What is detectable in experiments should not be the current jty itself because 
the QCD vacuum has fiuctuations of instantons and anti-instantons. In other words, 
the parity (P) and the charge-parity (CP) symmetries are broken only (spatially and 
temporally) locally at the topological excitation, but those symmetries are restored on 
average over fiuctuations. Thus, the CME-induced current jV is also a local object and 
the (ensemble or spatial) average makes it vanishing. In this way the CME is one of 
the candidates to give an account for the "Local Parity Violation" (LPV), if any, that 
might be observed in the relativistic heavy-ion collisions. 

The most relevant quantity to experimental data |160j is fiuctuation of the CME- 
induced current, that is, the electric-current susceptibility Xi [161] . The one-loop result 
in the small-frequency limit with the zero momentum limit first taken is 



•,2 



eB\ 



X. = ^ , (62) 

for QED loops (e^ should have been for quark loops), which does not depend on 
/is and comes from the Landau zero-mode alone, interestingly. There is an intuitive 
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argument to take a short-cut for the derivation of the above expression [161] . For this 
purpose let us consider the electric-current generation rate under strong E as well as 
-B, that is, 

d(eJy) _ ,, e\B\E 

which originates from the correspondence between the chirality generation and the 
particle production when fields are strong enough |162] . The same quantity can be 
expressed in the framework of the linear response theory as 



^ = - / d^xdV(^^^^jWx'))eA(x') 



dt J \ dt 

d^xd^x'e^jv{x)jv{x')) E , (64) 



where A{x) denotes a vector potential component parallel to B and Jy- From the first 
line to the second line above, we used E = doA. By equating these f l63|) and fl6^ . we 
can immediately find Xj given by equation (1621) . 



4-2. Implication to and from the QCD phase transitions 

Once the CME is confirmed in the heavy-ion collision experiment, it would signal chiral 
symmetry restoration. This is because, as we have seen in the previous subsection, the 
CME requires massless Dirac fermions and thus vanishingly small chiral condensate. 
It is not obvious what the CME can imply for deconfinement. Generally speaking, 
the sphaleron transition rate is proportional to by dimensional reason and real-time 
topological excitations become abundant at higher temperature |163[ I164| 1165] . In this 
argument, however, any feature inherent to deconfined gluons and quarks is not quite 
needed for the manifestation of the CME. In a solvable model in (1+1) dimensions, as 
we will discuss later in fact, the CME exists even though no deconfinement takes place. 

In the computation of the electric-current susceptibility there is a significant 
influence from the chiral phase transition. Actually, at finite /is, the divergent chiral 
susceptibility has a mixing with the current susceptibility Xj causing enhancement in 
Xj at the chiral phase transition |166j . It is not yet understood how the Polyakov loop 
dynamics and the deconfinement transition should affect Xj and any other observable 
sensitive to the CME. 

In principle the lattice-QCD simulation in a strong B field can clarify the effect 
of chiral restoration and deconfinement |123[ I124[ I167[ 1168] . There are also some 
effective model studies on the phase diagram modified by B and also the effects of 
the strong 6'-angle [I69l [ITOl [HIl HIS]. It is known by now that the PNJL and PQM 
model calculations with large B [173j are not consistent with the lattice-QCD results 
in which chiral restoration and deconfinement are locked together for any B [174j . This 
inconsistency could perhaps arise from the missing diagram as shown in figure [9] and the 
coupling to B in the Polyakov loop potential. Besides, the chiral model part may have 
a non-trivial dependence on the Polyakov loop through the fermionic interaction terms 
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[175[ 1176] . We also note that dynamical locking mechanism as discussed in [177] may 
play a role. At least, it should be quite robust that the chiral condensate is enhanced by 
B due to the magnetic catalysis |178l 1179] and thus the chiral transition temperature 
is pushed up accordingly. 

A careful comparison between results from the effective model approach and from 
the lattice-QCD simulation at large B should be very useful for the finite-density study 
on QCD matter. In the case at yUq > T it is difficult to impose any reliable constraint 
from the lattice-QCD simulation, while the finite-S simulation has no principle problem 
and one can check quantitatively if the effective model description is valid or inadequate. 
If it is insufficient to reproduce the lattice data at large -B, it is most unlikely that the 
same model can encompass the finite-density property of QCD matter either. 

One might have thought that the perturbative calculation of the Weiss 
potential ( fT2|) can be extended to the finite-S case. It is just straightforward to 
generalise the quark one-loop contribution ( !T8|) to the finite-S calculation. What is 
more difficult is the polarisation effect as in figure |9l This has been evaluated in the 
Lowest Landau-Level (LLL) approximation that will be explained in the next subsection 
for a quick derivation of the CME-induced current and susceptibility. In short summary, 
the polarisation effect on the Weiss potential has the following effect: The Weiss 
potential ( !T2|) comes from the integration with respect to two transverse gluons, which 
are physical degrees of freedom and unphysical longitudinal and ghost modes cancel out. 
In the LLL approximation one can show that only one of two transverse gluons acquires 
a screening mass proportional to a/ eB. The height of the potential, therefore, decreases 
up to a half of the original Weiss potential with increasing B. Qualitatively, a larger 
B tends to reduce the barrier at the confined state with $ = 0, strengthen confinement 
and delay the deconfinement phase transition. 

4-3. Dimensional reduction 

Under a strong magnetic field, in general, the transverse motion of charged particles is 
equivalent to the one in the harmonic oscillator. The energy level is then discrete due to 
the Landau quantisation. Spin-1/2 fermions have the Landau zero-mode which would 
dominate in the dynamics at energies below the scale ~ ^feB. Such a strong B enables 
us to use the LLL approximation and to drop the transverse motion completely. In this 
limit we can reduce the (3+l)-dimensional theory into a form of the (l+l)-dimensional 
one multiplied by the Landau level density. 

In Minkowskian space-time we use the metric g^^ = —g^^ = 1, g^^ = g^^ = 
and the 2x2 7-matrices which satisfy = 2g^'^. Chirality is characterised by 

^5 _ ^0^1 _ fiiag(^x, —1) in the chiral (Weyl) representation. Therefore the upper 
(lower) element of two-component spinor ip = [iPrjIPlY represents the right-handed 
(left-handed) particle. In (1+1) dimensions the particle-anti-particle difference is 
correlated with the chirality. That is, in momentum space, the right-handed component 
corresponds to a right-moving {p > 0) particle and a left-moving {p < 0) anti-particle. 
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One can understand the left-handed component in the same way, i.e. a left-moving 
(p < 0) particle and a right-moving (p > 0) anti-particle. 

In (1-1-1) dimensions the following relation among the 7-matrices plays an 
interesting role for the topological currents; 

where e^^ = —e^^ = —Eqi = £10 = 1, which relates the vector and the axial- vector 
currents. As usual, we can write the vector and the axial-vector currents as 

= i^ri^ , = ^7^'^ • (66) 

Using the relation (l65l) . we have a relation, = — e^'^jy, that is explicitly written as 

m 

Jv=jI Jl=f- (67) 

4-3.1. Topological currents in (1+1) dimensions: The relation between the vector and 
axial- vector currents is very useful because, as we will see here, it captures the essential 
feature of the CME-induced currents in (3-1-1) dimensions. 

Let us consider the anomaly relation in (H-l) dimensions. It is well-known that 
the axial anomaly leads to 

^mJs = ^^'"F.u = -E = -2qw , (68) 

where the electric field is E = F^^ in the standard convention. Note that there is 
no magnetic field but only the electric field E in (H-l) dimensions. We here defined 
the (l-l-l)-dimensional topological charge density as qw = —{e/2TT)E in accord to the 
convention. By integrating equation (1681) over space-time and assuming that the current 
falls sufficiently fast at spatial infinity, one can recover equation f l58l) easily. We can also 
prove that the topological charge, Qw = J (Pxqw{x), takes an integer number so that 
the boundary condition in the x-direction can be maintained. 

We also note that one can express equation fl68l) in the following form; 

d^j^ix) = -^(9^1 - 91^0) = -2d,K^ix) (69) 

with the (l-l-l)-dimensional Chern-Simons current density defined by K'^ = 
— {e/2Ti)e^^^ Ay. From this identification, the Chern-Simons number in this system is 
inferred as 



v{t) = J dxK%t,x) = ^ j dxA\t,x) 



(70) 



Combining these expressions with the relation jy = (where is the volume 
integral of jg), one can immediately write the vector current integrated over space as 

4 it) = N^it) = -2 I dt'dxqwit'.x) , (71) 
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assuming that was zero at the initial time {N^{t = — oo) = 0). This simple relation 
leads to the current at late time as given by 

4 = -2Qw . (72) 

This is nothing but the result expected when the spin is fully polarised in the (3+1)- 
dimensional CME at strong B (see equation ( !59l) ). Note that in (1+1) dimensions the 
spin is always fully polarised because there is only one spatial direction and thus the 
moving direction (either p > or p < 0) and the chirality of particles have one-to-one 
correspondence. Here equation ( 17Ti) physically means Ohm's law because the (1+1)- 
dimensional topological charge density is proportional to the electric field as seen in 
equation (|68l) . 

If the spatial component of the Chern-Simons current falls sufficiently fast, the 
topological charge is written as Qw = v{t = oo) — uit = — oo). Therefore, (the spatial 
average of) is the Chern-Simons number and the boundary condition of in the 
t-direction gives the topological winding number. Supposed that z/(t = — oo) = N^(t = 
— oo) = 0, the topologically induced current is written as 



4(t) = ~J dxA' 



(t,x). (73) 



If we identify —eA^ as the chemical potential /Xq (regarding the sign, remember the 
covariant derivative p° — eA^ and the dispersion relation p^ = Ep — for particles). 
Equation (165|) implies that eA^'j^ = eA^'y'^'y^ and thus —eA^ can be regarded as the 
axial (or chiral) chemical potential /^s. Therefore, one can reach a conclusion that 

4 = -[dxfi,, (74) 



TT _ 

which correctly recovers the (3+l)-dimensional CME-induced current fl60l) once this is 
multiplied by the Landau level density, eB/{2'K). That is, 

jv = — [in (1+1) dimensions] 

TT 

— J- jv = ^— — ^ ■ — , [in (3+1) dimensions] (75) 
zvr TT 

which coincides with equation fl60l) . 

Here, it is clear that the longitudinal gauge field A^, which is the Chern-Simons 
number in (1+1) dimensions, plays the role of the chiral chemical potential /is in (3+1) 
dimensions. We note, however, that there is an important difference; usually /is is 
introduced by hand as a constant, but in (1+1) dimensions A^ must have t-dependence 
to allow for nonzero Qw We can think of a concrete "instanton" configuration in (1+1) 
dimensions simply as 

M(t,x)='^L^-Et, (76) 

where we limit ourselves to the spatially homogeneous case and denote the spatial and 
temporal extents as L and T, respectively, and then we have 

Jy(^) = ^ = —t. (77) 



QCD matter in extreme environments 



32 



From this, again, if multiplied by the Landau-level degeneracy we can correctly recover 
the current generation rate given by equation fl63l) . i.e. 

d(ejy) e'^E 



dt TT 

d{e]v) \eB\ e^E 



[in (1+1) dimensions] 
[in (3+1) dimensions] (7J 



dt 271 71 

which coincides with equation ( l63i) . 

In the same way we can get a finite axial-vector current at finite quark chemical 
potential /iq. To see the anomalous nature in this case the important fact is that the 
relation between the density and the chemical potential is given by the quantum anomaly 
in (1+1) dimensions, that is, 

nq = -— , (79) 
vr 

which results from the (l + l)-dimensional anomaly. One can derive this expression 
directly from n = {ip'^{x)ip{x)) by inserting the gauge field gulator as 

limyO_j,2,.o ip'^iy) exp[— ie J dtA^]ilj{x). From this one can immediately find, 

J5 = y dxriq = y dxnq, (80) 

which represents the axial counterpart of the CME jl58j . This is again the anomaly 
relation exactly same as that in (3+1) dimensions once multiplied by the Landau level 
density eB/27i. 

4-3.2. Chiral Magnetic Effect in the Schwinger model: So far the arguments and the 
resulting expressions are quite general. From now on we shall go into the dynamical 
properties calculating microscopic quantities in a solvable (l+l)-dimensional model, 
i.e. the massless Schwinger model. The easiest way to accomplish a calculation in the 
Schwinger model is to use the mapping onto a free bosonic theory. In our case, however, 
the bosonisation rule is a bit more complicated than usual because we deal with not 
only fermionic fields (such as the chiral condensate) but also gauge fields (such as the 
electric field). So, the Lagrangian density of the corresponding theory should be 

lid'^9){d,9)-m.,id'^9){d,<P)-^( 
with the boson mass |181] . 

m^ = -. (82) 

71 

After integrating the 0- field out, we get a theory only in terms of the ^-field that is free 
(no interaction term) and has a mass m^. Such a scalar theory is usually used with the 
bosonisation rule |182] . 

= = -^e^^dj , (83) 

\/7l 



C = -{d^9){d,9) - m.,{d^9){d,4>) - -{d''4>)d\d,4>) (81) 



= ^7V^ = -^d^9 , (84) 

'7T 



ipifj = —cm^ : cos(2y7r6') : (85) 
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with the normal ordering : :. Now we remark that in the Lagrangian density fl8ip comes 
from the gauge field, = —e'^^'^dycj) (where (f) includes an instanton-like configuration 
~ \Et^ which does not satisfy the periodic boundary condition in the t-direction). Then 
the electric field takes a form E = d'^cj). Once we integrate the ^-field out from the theory, 
after the Gaussian integration in the functional formalism, equation (|8^ is replaced by 



777 p 

jt = — ^d^'e — ^d^(i) = --d^(t) . (86) 

^/7^ ^/7^ TT 



The anomaly relation is then derived as 

^mJs = --d^ = --E = -2qw , (87) 

TT TT 

which is fully consistent with the anomaly relation f l68p . 

In the same manner we can express the vector current in terms of to find, 

= -e^''d,<p = 2e^''%w . (88) 

The inverse Laplacian should be understood in frequency and momentum space (see 
equation (189|) ). It is easy to make sure that this result is fully consistent with the 
previous relation again. That is, after the spatial integration on (j) and qw in the above, 
the spatial derivative di drops off and the right-hand side simplifies as — 2/(9o for yU = 1 
component, that is just the t-integration. Therefore the right-hand side finally becomes 
—2Qw together with the spatial integration, and hence we obtain Jy = —2Qw- 

The above equation gives a microscopic structure of the current in more general 
cases with spatial modulation. In frequency and momentum space we can re-express 
this as follows; 

jy{uj, k) = qw{uj, k) . (89) 

This is an interesting relation. If a; — is taken first, we see that jy(0, k) is vanishing. 
To get the CME-induced current and the non-zero chiral magnetic conductivity, it is 
necessary to take the zero-momentum limit in the order of /c — )■ first and then w — )■ 
later. This observation is in fact consistent with the result of the one-loop calculation 
of the conductivity [147] . 

We point out that the structure of equation ( l89il naturally appears from the 
transverse projection. That is, after the one-loop integration with the gauge potential 
source in momentum space, the well-known result reads; 

j^iu, k) = - (^g^'' - ^A,{u, k) (90) 

with q = {u, k) and = u;^ — A;^, from which one can easily find that 

3v{^,k) = —^^-A\u:,k). (91) 

Because qw = {e/T^)d'^A^ , one can substitute = i{2'K / e)qw /oj for A^ above, and one 
can then check explicitly that the above expression is equivalent to equation 
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From the equivalence to the bosonised theory it is very easy to read the electric 
current-current fluctuation too. To this end one should integrate the 0-fleld flrst, and 
then what remains is a free massive scalar theory in terms of the ^-fleld alone. Then we 
trivially get, 

-y) = e'{j\x)j\y)) = ^o^^o^ {e{x)e{y)) , (92) 
or in momentum space one can express this as 



^2 _ f^2 _ ^2 _|_ 



(93) 



At a flrst glance this expression looks different from equation (15^ . This is because the 
above expression f l93|) is a result after resummation of the bubble-type diagrams, while 
equation fl62l) is the result of the one-loop order. Roughly speaking, appears in the 
denominator of equation ( l93l) as a result of inflnite insertion of the polarisation diagram. 
This indicates that one can extract the one-loop result from the leading-order Taylor 
expansion of equation ( l93l) in terms of m^. Such a procedure actually leads to 

^2 



Therefore, 



one— loop , 
Xj I 



one— loop 



u, k) 



A;2 



— )■ 



e 
vr 



— (at A; ^ 0) . 



(94) 



e 

TT 



one— loop 



\eB\ 



e 

TT 



[in (1+1) dimensions] 
[in (3+1) dimensions] 



(95) 



which again coincides with the previous result (1621) . 

In this way the Schwinger model is of great use to understand the topological 
properties probed by the magnetic fleld. In reference [183 ] microscopic calculations 
to demonstrate how the dimensional reduction occurs in a way consistent with the 
momentum conservation are given with a more detailed result for the polarisation tensor 
in (1+1) dimensions embedded in (3+l)-dimensional gauge flelds. 



5. Summary and Outlook 

In this review some of the theoretical approaches to QCD matter in extreme 
environments have been picked up. This direction of physics is strongly motivated by 
relativistic heavy-ion collision experiments. Furthermore, extreme environments such 
as the high temperature, the high baryon density and the strong magnetic fleld would 
enable theorists to attack QCD problems in a treatable way. 

Theoretical and experimental researches on flnite-T QCD have achieved the level 
of the precision science, whereas the flnite-density study of QCD is still controversial. 
Theoretical approaches cannot escape from huge uncertainties, and only the forthcoming 
experimental data will be able to impose constraints on many possibilities proposed from 
the theoretical side. 
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New physics opportunities provided by the strong magnetic field created in the 
heavy-ion colhsion are quite intriguing. A deeper understanding in this direction would 
be helpful for the finite-density study as well. This is because gluons can be coupled to 
both the magnetic field effects and the density effects only through the quark polarisation 
processes. 

There are many interesting subjects that we had to miss in this article. Let us 
look quickly over some of them. We did not discuss the recent developments in the 
functional approaches to the QCD phase diagram based on the renormalisation group 
flow |184j and the Dyson- Schwinger equation [185[ 1186] . This approach is the most 
promising among others to attack the problem of the QCD phase diagram from the 
flrst-principle technique. It has been understood how quark conflnement is realised in 
terms of the gluon and ghost propagators |187j . Future extensions to the three-flavour 
case without uncertainty that stems from the (T-dependent) strength of the U(1)a 
anomaly would be desirable if possible. The inhomogeneous chiral condensate (chiral 
density wave) is the key concept that may reconcile various states of matter in the 
baryon-rich regime; the QCD critical point, the QCD triple point |188j . Quarkyonic 
Matter and even the dimensionally-reduced state in strong magnetic flelds |180] . The 
relation between the chiral density-wave state and the Polyakov loop dynamics would 
be a challenging problem too. 

Intense magnetic flelds have opened a new direction of physics in the heavy-ion 
collision. The situation realised as a result of the LLL approximation is similar to 
that near the Fermi surface at high density, i.e. pseudo-(l+l) dimensionality causes 
peculiar phenomena such as superconductivity |189j and sound modes |190[ 1191] that 
may have a connection to the so-called Tomonaga-Luttinger liquid in which no quasi- 
particle excitations but only sound modes exist. 

All these developments and new possibilities are waiting for further investigations. 
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